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We study models with a generalized inhomogeneous equation of state fluids, in the context of 
singular inflation, focusing to so-called Type IV singular evolution. In the simplest case, this cosmo¬ 
logical fluid is described by an equation of state with constant w, and therefore a direct modification 
of this constant w fluid, is achieved by using a generalized form of an equation of state. We inves¬ 
tigate from which models with generalized phenomenological equation of state, a Type IV singular 
inflation can be generated and what the phenomenological implications of this singularity would 
be. We support our results with illustrative examples and we also study the impact of the Type 
IV singularities on the slow-roll parameters and on the observational inflationary indices, showing 
the consistency with Planck mission results. The unification of singular inflation with singular dark 
energy era for specific generalized fluids is also proposed. 
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I. INTRODUCTION 

One of the most astonishing surprises that the scientific community experienced in the late 90’s was the observa- 
tionally verified late-time acceleration of the Universe [l|, an observation which was contrary to any up to that date 
perception or expectation for the Universe’s late-time evolution. Since then, considerable amount of research was 
conducted towards the understanding and explanation of this late-time acceleration. In addition, a desirable feature 
that a complete theory of cosmological evolution should have is the description of the early-time acceleration, known 
as inflationary era ii, and of the late-time acceleration using the same theoretical framework. One of the most 
successful approaches that describe late-time and early-time acceleration in unified manner Q is provided by F{R) 
theories which serve as modifications of the standard Einstein-Hilbert gravity. It is very interesting that modi¬ 
fication of gravity may be considered as addition of (geometric) terms to the phenomenological equation of state of 
corresponding generalized fluid. In fact, it is quite well-known that generalized (imperfect) fluid EoS Q which contain 
inhomogeneous terms may successfully describe the universe acceleration. The homogeneous modifications assume 
the inclusion of terms in the EoS which depend on the effective energy density, and the inhomogeneous modifications 
assume terms that depend explicitly on the Hubble rate or its higher derivatives, and as was demonstrated in Q, in the 
context of inhomogeneous phenomenological EoS theories, it is possible to describe even phantom evolution without 
introducing a scalar field with negative kinetic energy. Apart from this appealing feature of phenomenological EoS 
theories, a strong motivation to use and study these comes from the fact that inhomogeneous terms may be understood 
as the time-dependent bulk or shear viscosity Q and also symmetry considerations indicate such modifications of EoS 
3 . Such generalized fluids may be used for the construction of inflationary era (see, for instance [I3|). Einally and 
most importantly, it can be proven that in the context of modified gravity theories, the resulting EoS of the geometric 
fluid is modified in the fashion we described above 3 111 • 

It is common in the context of these theories, however, that various types of finite time singularities frequently occur. 
It is remarkable that such finite-time singularities are universal, as they may occur after early-time acceleration as 
well as after late-time acceleration. The purpose of this article is to study the Universe’s cosmological evolution using 
a phenomenological EoS for the dark fluid, that leads to, or is responsible for, a Type IV finite time singularities. The 
interest to this specific singularity is owing to the fact that the Universe’s evolution continues smoothly after passing 
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through this singularity! This argument is strongly supported by the fact that no geodesics incompleteness occurs for 
sudden and this types of singularities, so in principle the evolution is not abruptly interrupted by the singular behavior 
of some physical quantities, such as the higher derivatives of the Hubble rate. The study of realistic cosmological 
singularities was initiated sometime ago [HI- The very interestin g fin ite time future singularities are known as sudden 
singularities or Type II ones (for later studied on these, see |12 h 141 |'). Also the finite time cosmological singularities 
were firstly classified in , with the classification taking into account whether physical quantities like the scale factor, 
the effective energy density, or pressure, are finite or infinite at certain spacetime points. The Type IV singularity is 
the most “mild” singularity among the four types of singularities, with the mild referring to the fact that these are 
not crushing type singularities. For recent studies on these singularities see [l6l - [T8l| . where details can be found for 
the various types of cosmological evolution that these can generate or modify. 

The purpose of this letter is to demonstrate at first time that singular inflation induced by generalized fluids is quite 
possible. The unification of singular inflation with singular dark energy era in the universe filled with generalized 
fluid is also proposed. The paper outline is: In section II, we present some essential information for the theoretical 
framework we shall use and we study in detail the how a singular evolution can be produced by generalized EoS 
models and we also provide some concrete examples. In section III, we address the slow-roll behavior issue of these 
generalized EoS theories and in section IV the conclusions follow. 


II. THEORETICAL FRAMEWORK AND ANALYSIS OF SINGULAR GENERALIZED EOS MODELS 

To start with, consider a perfect fluid coupled with the standard Einstein-Hilbert gravity, with the corresponding 
Friedmann-Robertson-Walker (FRW) equations being equal to, 

p = + , ( 1 ) 

where it is assumed that the background metric is a flat FRW metric of the following form, 

ds^ =-dt^+a^(^)^dx^ (2) 

i 

In Eq. ([T|), H denotes the usual the Hubble rate H{t) = d{t)/a{t). For a given cosmological evolution in terms of the 
Hubble rate H — H{t), the right hand side of the two equations appearing in Eq. ([T]) are solely functions of t, 

P = fpit ), P = fpit) ■ (3) 

Then, by solving the first equation with respect to the cosmic time t, we get, t = fp^ (p), and by substituting the 
resulting expression into the second equation of Eq. dSD, we obtain the functional form of the EoS, namely, 

P = fp ifp^ (P)) • (4) 

In this way we can obtain the exact form of the effective EoS, which we denote as Weff ■ It is worth recalling in brief the 
essential features of the homogeneous and inhomogeneous phenomenological EoS theories, the complete presentation 
of which can be found in Q . In the context of inhomogeneous phenomenological EoS theories the EoS is modified as 
follows. 


P = -p-f{p) + G{H), (5) 

and when the function G{H) is zero, this corresponds to the homogeneous EoS theory. A much more general functional 
dependence of the effective pressure as a function of the effective energy density and the Hubble rate, is given by the 
following functional form, 


P = f{p,H). (6) 

In Eqs. dU and ([6]), the functions / and G are arbitrary functions of their arguments. 

The focus in this letter is to study the phenomenological EoS theories of Eq. m, for a Type IV cosmological 
evolution as it occurs in the inflationary epoch. It is worth recalling in brief the classification of finite time cosmological 
singularities, as in Ref. [l5l| : 

• Type I (“Big Rip”) : When t —>• 4, the scale factor a, the effective energy density pes, and the effective pressure 
PeS diverge, namely a —>■ oo, Peff oo, and |peff| oo. For details on the Big Rip finite singularity, the reader 
is referred to [l^ [l^ [2^ . 
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• Type II (“sudden”): When t ^ tg, although both of the scale factor and the effective energy density are finite, 
that is, a —>■ Og, Peff Ps, the effective pressure diverges, namely |peff| oo. See [3 [11 for further analysis 
on this type of singularities. 

• Type III: When t —>• tg, although the scale factor is finite, a —>■ Og, both of the effective energy density and the 
effective pressure diverge, peff —t oo, |peff| —t oo. 


• Type IV: When t —>■ tg, all of the scale factor, the effective energy density, and the effective pressure are finite, 
that is, a — >■ Og, pes — t Ps, IPefrJ— t Ps, but the higher derivatives of the Hubble rate diverge. For a detailed 
account on this singularity, see [ll. 

Here the effective energy density peS, and the effective pressure Peff are defined by peff = (3 /k^) H^, and Peff = 
— (1 /k^) A quite general form of the Hubble rate that can lead to Type H and Type IV singular 

cosmological evolution, is the following. 


H{t) = fi{t) + / 2 (t) {ts - , 


(7) 


where the functions fi{t) and / 2 (t) are considered to be smooth and differentiable functions of t. The Type H 
singularity occurs when the constant parameter a is restricted to take the values, 0 < a < 1 , while when a > 1 , 
the cosmological evolution develops a Type IV singularity. In the following, and in order to avoid inconsistencies, we 
assume that the parameter a appearing in Eq. o, takes the following form, 


a = 


n 

2m + 1 


( 8 ) 


where n and m are arbitrarily chosen positive integers. In addition, a modified version of the Hubble rate appearing 
in Eq. 0 , is given below. 


H{t) = h{t)+h{t)\tg-tr, ( 9 ) 

in which case, a can take more general values, without assuming the restricted form given in Eq. (jS]), but we study 
here only the case for which a has the form given in Eq. ([ 8 ]) . In the following we shall investigate which homogeneous 
and inhomogeneous EoS can generate the cosmological evolution appearing in Eq. o, with special emphasis given 
in the Type IV singularity. In general, it is difficult to find an explicit form of the EoS, for complicated forms of the 
functions fi{t) and f 2 (t), so let us start with a simple example describing a Type IV evolution, with the functions 
fi{t) and / 2 (t) being chosen as, /i(t) = 0 and f^it) = /o, where /o is an arbitrary positive constant. By using these 
values and by substituting in Eq. o the Hubble rate of Eq. o, the effective energy density p and the effective 
pressure p read. 


p=^{tg- , p = —\ {ts - tf“ + 2afo (ts - t)“ , (10) 

By using the first equation of Eq. (HI, we can solve it in terms of ts — t and therefore by substituting in the effective 
pressure expression, we obtain in explicit form of the EoS, which reads (see also [3), 

p = -p-2-3-^aK-^f^-^p^ . ( 11 ) 

The exact type of the finite time singularities that may appear in this cosmological evolution may be found by using 
the explicit forms of the effective energy density p and of the effective pressure p, as these appear in Eq. (nni). By 
introducing the parameter a = 3 ^, the cosmological evolution of the phenomenological EoS (ITT]) , as described by 
the energy density and pressure appearing in Eq. (O, has the following singularity structure, 

• In the case that a > 1, or in terms of a, —I < a < —1, a Type HI singularity occurs. 

• In the case that 1 < a < I, that is, a < —I, a Type I singularity appears. 

• In the case that 0 < d < 5 , or equivalently, a > 1, a Type IV singularity occurs. 

• In the case that d < 0, or in terms of a, —^ < a < 1, a Type H singularity occurs. 
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The form of the EoS appearing in Eq. (ED, can be viewed as a homogeneous phenomenological EoS theory of the 
form dS]), with f{p) being equal to, 


/(p) = -2.3-^a«-^/o^-V^. (12) 

and G{H) = 0. In addition, the EoS appearing in Eq. dill) can also be viewed as an inhomogeneous phenomenological 
EoS theory of the form (O , by using the fact that p can be written in terms of the Hubble rate H'^ in the way dictated 
by Eq. (HD. By doing so, the EoS can be written in the following form, 


j* 2 —a 

P=-p-^fo 


H- 


which is of the form given in Eq. dSD, with f{p) = 0 and with G{H) being equal to, 

G{H) = . 


(13) 


(14) 


We may also consider the following model, 


H{t) = ho 



(15) 


where /iq, fo, H, a are constants and we assume ho > 0, n > 0 and also that 0 < a < 1. a > 1. When t —>■ ±oo, 
we find H{t) becomes a constant H{t) —>■ ho and when t ^ to, a. Type IV singularity occurs, since the Hubble rate 

behaves as, H{t) ^ ho ■ The corresponding effective energy density and effective pressure for the Hubble rate 

(fTSI) are equal to, 



P = - 







Consequently, by making use of Eq. o, we obtain the following EoS, 

__ 2aho / K^p\ / / K^p\ 

^ ^ VS/iq/ \v3/io/ 



Now consider the case in which the Hubble rate is given by, 

H{t) = fo{t - ti)“ + co{t - 12^ , 


(16) 


(17) 


(18) 


with Co and /o constant and positive parameters, and a, /3 > 1, so that the cosmological evolution has two Type IV 
singularities at f = and t = t 2 . We may choose ti to be at the end of the inflationary era and t 2 to be at late-time. 
In principle, it is quite difficult to obtain the exact form of the EoS for the Hubble rate (fTSl) . since it is quite difficult 
to solve the equation p ^ H‘^(t) with respect to t. However we can find an approximate form of the EoS, near the two 
Type IV singularities. Before going into the details of this approximation, we quote here the cosmic time dependence 
of the effective energy and of the effective pressure for the Hubble rate (fT51) . which are. 


_3/o(t-ti)^“ Qcofo{t-tir{t-t2)^ 3c^{t-t2ff^ 

P o' 9 ' 0 ^ 

_ 3/|f(<-<i)2“ 6co/o(t - ti)“(t - ^ 2 )^ ‘icl{t-t 2 y^ 2/o(t - fi)“^+“a 2co(t - t2)“^+^/3 

P ~ ^2 ^2 „2 ^2 ^2 

rp Ay A; A< A< 

So at the vicinity of the early-time Type IV singularity, these read, 

_3co(t-t2)^^ _ 3co(t-t2)^^ 2coit -12)~^~'~^ 13 

P ~ ..2 >7* .9 ,.9 5 


(19) 


( 20 ) 
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so the EoS takes the approximate form, 


P = -P- 


2co/3 


Wo) 


( 21 ) 


The resulting physical situation is quite appealing, since it seems that the late-time singularity controls the early-time 
EoS for the inhomogeneous phenomenological EoS theory that generates the Hubble rate (TT^ . near the of course 
early-time singularity. In fact, the resulting inhomogeneous phenomenological EoS near the early-time singularity is 
controlled solely from the late-time singularity, since the terms ^ (t — ti) vanish for a Type IV early-time singularity. 
Of course, for other types of singularities, this may not occur, so only the Type IV case has this interesting feature. 
The same applies for the late-time singularity, so near the late-time Type IV singularity the effective energy density 
and effective pressure are, 


P = 


icl{t - hf 


P = 


3cg(t-ti)2“ 2co(t-ti) 


— l+a 


SO that the corresponding EoS reads. 


2coQ! 

p= -p - — 



( 22 ) 


(23) 


Therefore, the resulting picture is that the EoS near the late-time singularity is solely controlled from the early-time 
Type IV singularity, which is a quite interesting feature. So by suitably choosing the parameters, this effect can be 
quite large or even negligible. For a relevant study in which this phenomenology also occurs, see also (l^ . 

Another interesting model with the property that it provides a unified description of inflation at early-time with a 
late-time acceleration evolution with a Type IV singularity occurring at late-time. The Hubble rate of this model is 
given below. 


H{t) = 


h 


+ 


h{-t + t2Y ■ 


(24) 


Notice that the late-time singularity can even occur at present time. For a detailed account on the consequences of 
this Type IV singularity occurring at present time, see [1^. In Eq. (IMl) . the parameters a, /3, fi, / 2 , and /a are 
chosen to be positive constants so that we ensure that H{t) > 0. Then, the effective energy density p and the effective 
pressure p are equal to. 


fi , + 




t^ + 4 


+ /3(—i \ , 


P = - 


fi , + 


\/t‘^ + t'l 


t^ + 4 


+ f3i-t + t2)^ 


fit 


^f2tH-t + tir , 2f2t{-t + tiy 


I {t^ + tlf 


+ d 


(25) 


The cosmic dark fluid scenario with effective pressure and energy density as in Eq. (1^51) can describe a plethora of 
cosmological evolutions. For example, ti may correspond to early-time and t 2 at late-time, so if a > 1 and /3 > 1, a 
Type IV singularity occurs at both early- and late-time, as in the previous example. Let us briefly repeat the same 
analysis as in the previous case, in order to find an analytic approximation for the EoS near the singularities. Consider 
first the case that t ~ ti, so the physical system is considered to be near the early-time Type IV singularity. Then, 
the effective energy density and pressure become approximately equal to. 


3/i &hh{-t + t2y + 

2/if _ ifl _ 6/i/3(-t + t2)^ _ 3/|(-t + t2)^^ 2/3(-t + f2)-^+^/3 

^ ( t 2 + i 2 ) 3 / 2^2 ( t 2 + t 2)^2 ^^2 ^ 2^2 ^2 ^2 


( 26 ) 
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and as is obvious, the late-time Type IV singularity controls the effective energy density and effective pressure of the 
cosmological dark fluid, near the early-time Type IV singularity. In order to see how the late-time singularity controls 
the EoS of the phenomenological theory, let us further simplify the expressions appearing in Eq. (1261) . by taking into 
account that t <^t 2 , in which case the effective energy density and effective pressure read. 


p ~ 


3/|(-< + kff 


3/|(f 


P~-- 






(27) 


so that the EoS reads, p ~ —p. Also notice that, owing to the appearance of the parameter t 2 , the late-time singularity 
drastically affects the early-time EoS and acceleration. The latter, is due to the fact that the early-time evolution is 
nearly de Sitter acceleration (since p ~ —p)- 

Conversely, at late-time and near the future Type IV singularity, that is, when t ~ t 2 , the corresponding effective 
energy density and pressure of the phenomenological EoS theory are. 


3/f 6/o/it"(-t + ti)“ + 

^ (t2-pt2)K2 (^4 + ^4) ^2 (i4_^i4)2^2 

2/it _ 3/f 8/ot5(-t + ti)" _ 4/ot(-t + fi)" _ + 

^ “ (i2 ^2)3/2 ^2 (t2+t2)^2 (^4 ^4)2 ^2 (^4 + ^4) ^2 ^^2 + ^2 (^4 + ^4) ^2 

3/gt4(-t + fi)^“ 2/ot"(-t + ti)-i+"a 

(t4 -p t4)2 ^2 (t4 + tl) k2 


SO practically, we could say that the early-time Type IV singularity controls the late-time behavior of the phenomeno¬ 
logical EoS theory. Let us see however the extent of this control explicitly. Since we assumed that t ^ ti, the effective 
energy density and pressure appearing in Eq. (12811 are very much simplified, and by keeping only leading order terms, 
these become approximately equal to. 


(t4+t4)2^2 ’ P (i4 + i4)2^2 


(29) 


and hence it easily follows that in this case, the EoS is p ~ —p, so we have late-time de Sitter acceleration with 
iCefi — —1. Therefore, the contribution of the early-time singularity is not so important in this case. Another 
interesting example with two different Type IV singularities is described by the following cosmological evolution. 


H{t) =/o + c(t-ti)“ (t-t 2 f , 


(30) 


in which case, if a,/3 > 1, Type IV singularities occur at t = and at t = t 2 - The corresponding effective energy 
density and pressure are equal to. 


P = 

p = 


2 

2 


3c2 

6c/o(-t -b tl)“(-t + t 2 )^ 

3/o"(-i + ii)"“(-i + l2)"'^ 




2/o(- 

-t + + t2)^a 1 

2/o(~l + li)“(~^ + I 2 ) 


At both singular points, the effective energy and pressure are exactly equal to, 

3c^ 3c^ 

P y 7 P y : 


(31) 


(32) 


which implies that the EoS is exactly equal to minus one, since p = —p. We have to note that in all the above 
paradigms, the results can drastically change, if the singularities we considered are not Type IV. 

The converse procedure is possible to produce a Type IV singular cosmological evolution, from a given phenomeno¬ 
logical EoS theory. Indeed, consider a homogeneous phenomenological EoS theory, with the EoS being of the form. 


P = -P + f{p), 


(33) 


with f{p) being equal to, 


f{p) = Ap^ . 


(34) 
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Then, by using the energy conservation law, we easily obtain that. 


P= 



Correspondingly, the effective pressure as a function of time reads, 

P = - 

X 

so when 0 < a < i, a Type IV singularity occurs in the cosmological evolution, at t = to- 


-2a 


(^) - to) {An) 1-2. (^(3) _ ig) 1-^2. (Ak) 

(1 — 2a)‘^K^ 

((!)— ^ ((|)~ ) 

(1 - 2a)2K2 


(35) 


(36) 


III. SLOW-ROLL PARAMETERS FOR SINGULAR GENERALIZED EOS MODELS 


In all the above cases, we investigated the general behavior of various homogeneous or inhomogeneous phenomeno¬ 
logical EoS theories, in the presence of a Type IV singularity during the cosmological evolution. However, as was 
stressed in [l^ , the presence of a Type IV singularity during the inflationary era may have dramatic consequences on 
the inflationary observational indices. It is therefore compelling to study the behavior of these observational indices 
in the presence of a Type IV singularity. We shall stay at a qualitative level analysis, but we shall attempt to compare 
some results with the recent Planck observational data [2l|. This analysis however is strongly model dependent, so 
in this letter we will just highlight the most important qualitative implications of the Type IV singularities. 

The full analysis on the derivation of the slow-roll parameters for a perfect fluid phenomenological EoS theory can 
be found in Refs. 0, and the key point is that the FRW equations appearing in Eq. m can be expressed in terms of 
the e-folding number N, as follows. 




p{N) + p{N) = 


2H{N)H'{N) 


(37) 


where H'{N) = dH/dN. Assuming that the EoS has the following general form, 


p{N) = -pmatter(A^) + f{p[N)) , 
the second equation in (l37l) . takes the following form, 

The energy density and the pressure satisfy the conservation law, 

p'{N) + 3H{N) {p{N) + p{N)) = 0 , 


(38) 


(39) 


(40) 


where p'{N) = df{p{N))/dN and in view of Eq. (IMl) . this conservation law becomes, 


p'{N)+3f{p{N)) = 0. 


(41) 


Combining Eqs. dnj and (1551) . we obtain, 

4 [{H'{N))^ + H{N) + H"{N)] = 3f{p)f{p) , (42) 

where in this case f'{p{N) stands for f'{p{N) = df{p)/dp. Consequently, we may express the slow-roll parameters 
as functions only of p{N) and f{p{N)) or equivalently as functions of H{N) and it’s derivatives (for details see @). 
For the present purposes we shall use the slow-roll parameters as functions of H{N), but later on we shall use the 







explicit form of the function f{p(N)), in order to have comparison with observational data coming from Planck 
The slow-roll parameters e, rj and ^ are given in terms of H{N) as follows, 


e = — 


H{N) 

AH'{N) 


6H'{N) 

H{N) 


H"(N) 

HW 


( H'(N) y 

yrnw)) 


1 2 


1 / H'{N)\ 

^-"2 


H'{N) 

H{N) 


9H'{N) 3H"{N) 1 fH'{N)\ 


H{N) H{N) 


+ 


V HiN) ) 


e = 


GH'jN) H"iN) 
H{N) H{N) 


f H'(N) y 
V H(N) ) 


H{N)H''{Ny 5H"'{N) 


3H{N)H'"{N) 9H'{N) 


1 l'H"{N)y 3H''{N) H"'{N) 

~ 2 \ H'{N) ) H'{N) H'{N) 

2H(N)H"(N)H"'{N) AH''{N) 


+ 


H'{NY 


H'{N) 


H'{NY 

3H{N)H"{NY 

H'{Nf 


H{N) 

'H"{N) 
H'{N) 


+ 


h'{ny 

15H"{N) 

H'{N) 


H{N) 

H{N)H""{N) 

H'{NY 


(43) 


It is convenient to have the slow-roll parameters as explicit functions of the cosmic time, in order to examine the 
various forms of the Hubble rate we presented in this letter, that actually lead to a Type IV singularity. In this 
way, we will explicitly study how the singularity that occurs at a finite time tg, affects the slow-roll parameters and 
therefore the approximation itself. These expressions are given below. 



r = 


i/4 

H'^H 

i/2 

3// 

// 

2i/4 

' i/s 

2//2 

HH 

i/2i/ 


2H 1 

4//2 

HH 

3//2 

ii2 

H^H 

i/2i/2 

i///3 

i/3 


-k 


#3 


HH^ nm I 


+ 


(44) 


Assume for the moment that the Hubble rate H (t) is given by the very general form of Eq. ([71) , and therefore when 
a > 1, a Type IV singularity is realized. Also assume that the function fi{t) is constrained in such a way so that 
fiYs), fiYs), and flits) do not vanish. Then, the slow-roll parameters at the vicinity of the Type IV singularity 
t ^ ts, behave as follows. 


hits 


4/i(ts 

fljts) 

4/i(i. 

-h (3 + 


flit.)- + TTItlFj + 


flits) 

F 


-2 


/ 2 (tg)a(a - 1 ) its - t)“ ^ (3 -k 


when a > 2 
when 2 > a > 1 


flits) \ 

2 V" tfffJ 

/ Gfljts) J_ fljts) 


V hits) 


2 / 1 ( 4 .)“ 
-1 


flits) 

/TIFF 


flits)' 

2/l(t.)“ 


/l(4d /l(4.) 

hits)^ 


/l(4d 

2/1 (t.)" 


3 flits) 


flits) 


flits)flits) flits)^flits] 


1 / 6 / 1 ( 4 .) , fljts) 

4 VbViF /i(4, 
9 / 1 ( 4 .) 


-1 


1 /q I fljts) \ “ /2a(a-l)(a-2) /, _ ,i 

2 + fiit,)^ ) /i(4.) = /i(4.) 

fljts) \ 

hitiv) 

i . 


a-3 


when a > 3 
when 3 > a > 1 


,/l(4.)/l(4. 

3 / 1 ( 4 .)'' 

■ AdlT' 


1 / V 7 IW, 

37 1(4.) , 27,(4.) 4/1(4.)“ 

fljts)'" /l(4.) = /l(4.) /i(4.)2/i(t.)“ 


fljts)' , f 1 ( 4 .) 

/l(4.)/l(4.)* /l(4_.^/l(4.)“ 


) 


hits) f 1 ( 4 .) 

hits) hits)" 


1 / 6 fljts) I fljts) \ fo_i_ fljts) A /2(4.)a(a-l)(a-2)(a-3) /, _ ,i 

4 [fljts)^ + fljthy flit.)-) flits)flits)'" 

ifo , fiits)\ ^ f2its)^a-{a-l)-(a-2)(a-3) f, 

+ /i(t.)4/ip,)^ 


a—4 


\2a-6 


when a > 4 
when 4 > a > 2 
when 2 > a > 1 


(45) 


If flit) is assumed to be smooth, the slow-roll parameter e blows up when 2 > a > 1, which means that it develops 
a singularity, while it is regular for a > 2. In addition, the slow-roll parameter rj, blows up for 3 > a > 1, while 
blows up in two different ways when 2 > a > 1 and when 4 > a > 2. Therefore, when a > 4, the slow-roll parameters 
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contain no singularities at least in the vicinity of the Type IV singularity. It is worth presenting some illustrative 
examples in order to further scrutinize the behavior of the slow-roll parameters. In the case that H{t) = foit — 
the slow-roll parameter e, becomes equal to, 

_ fojt - + 6t- 6ts + g)^ , . 

4(3/o(t-4)i+“+af 


while the slow-roll parameter rj is equal to, 

{t — (—-I- 8tts — 4:t1 + At'^a — Sttsa + -\- 2ttsa'^) 

^ 4/o (3/o(t — -I-a) 

(t — + 2o? — 2a^ — 6fo(t — 1 -I- 30;) -I- /q (t — + 2(—I -I- a)a)) 

4/o (3fo(t - -F a) 

Finally, the slow-roll parameter as a function of the cosmic time t, is equal to. 


^2_{t-is) ® ^“(-I-I-a)(-l-I-6t - 6ts-I-a) 

^ " 4/oM3/o(t-t.)i+“ + a) 

X (5(t - tsf{-l + aY + 3/o (f - tsf°‘[-2 + aYi-l + a)a + 3/o(t - ts)^+“(l -I- 2 q;)) . (48) 

By looking Eqs. and (H51) . we can easily see that the slow-roll parameters have singularities at the point 

t = ts, where the Type IV singularity occurs, a fact that we also stressed in the general example we presented earlier. 
The singularity in the slow-roll parameters can be viewed as rather unwanted features of the theory, or these can 
indicate a strong instability of the dynamical system that describes the cosmological evolution. Work is in progress 
towards the latter possibility. 

As a hnal task we shall try to compare the results of our analysis with the recent Planck data [2l|, by suitably 
choosing the phenomenological EoS function. For the purposes of our analysis, we shall express the observational 
indices as functions of the e-folding number N. We shall consider three important observational indices, namely the 
spectral index of primordial curvature perturbations ris, the scalar-to-tensor ration r and the running of the spectral 
index Os, which as was evinced in Ref. these are written in terms of the e-folding number N as follows. 


n, - 1 = - 9p(7V)/(p(7V)) 


( f'{p{N))-2 y 6p(V) i f{p{N)) 

\2p{N)-f{p{N))) 2p{N)-f{p{N))\ p{N) 


, 1 5f{piN))f{p{N)) , ff{p)\ 

+ +f{piN))-- -- +[^)) 


r =2Ap{N)f{p{N)) 
a. =p{N)f{p{N)) ( 


+ f{p{N))f"{p{N)) - 2 


~np{N))-2 y 

2p{N) - hp{N))) ’ 
f{p{N))-2 Y 


1 P'{N) 

\piN)J ' 3/(p(iV)) 

f{p{N))f{p{N)) ( f{p{N)) \ 


piN) 


-54p(iV)/(p(iV)) 


2piN) - fipiN))) 


72p{N) 


2p[N) - fipiN)) 




s 


2piN)-fipiN))) ripiN))-2 


J 2 


\ P{N) ) 


(49) 

(50) 


(51) 


where the detailed functional form of Ji and J 2 is given in the appendix. In order to have a qualitative idea of the 
behavior of the observational indices in terms of the phenomenological EoS of Eq. dSSD, we shall study a not very 
sophisticated model, but simple nevertheless, that leads to a Type IV singularity. Suppose that the EoS is given by 
Eq. (1331) . with fip) = Ap°‘ and in order to proceed we have to express the function /(p), in terms of the e-folding 
number N. This can be easily done, since the scale factor in this case is equal to, 

ait) = , 


(52) 
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so the effective energy density as a function of the e-folding number N is given by, 

p = - a)A)^ . (53) 

Recall that the function f{p) can generate a Type IV singular evolution when 0 < a < as we demonstrated below 
Eq. (1551) and therefore, the fraction f{p)/p can be chosen to satisfy the constraint. 


— < 1. (54) 

P 

In this case, the observational indices can be very much simplified, and can be approximated by the following expres¬ 
sions 0, 


Us — I — 6 


f{p) 

p(7V)’ 


r ~ 24 


f{p) 

p{N)^ 


Combining Eqs. (1531) and (1551) . with the constraint (l54ll holding 
observational indices read, 


true, the resulting approximate expressions of the 


2 ^,8 ^ 1 
^ “ V(1 - a) ’ V(l-a)’ “ “V2(l-a)2 ’ 

The 2015 Planck report [2l|, restricts the values of the observational indices as follows, 

= 0.9644 ±0.0049, r<0.10, =-0.0057 ± 0.0071, 


(56) 


(57) 


which can be relaxed if someone assumes a scale dependence of the scalar and tensor spectral indices. By using the 
values {N,a) = (60,1/20), the observational indices become equal to, 

71^- 0.96491, r- 0.1403, =-0.000307. (58) 


Therefore, there is concordance with the spectral index of primordial curvature perturbations, but the scalar-to-tensor 
ratio and the associated running of the spectral index constraints are not satisfied. However, in principle concordance 
can be achieved, if a more sophisticated model is considered, instead of the simple model we used here, just for 
expositional purposes. 


IV. CONCLUSIONS 

In this letter we investigated various phenomenological EoS models that can generate a Type IV singular cosmo¬ 
logical evolution. The motivation for considering phenomenological EoS models is coming from the fact that the 
Universe’s EoS seems to be non-constant and it may slightly cross the phantom divide at late-time. In addition, the 
Type IV singularity is the mildest possibility of singular evolutions, since it is of non-crushing type, and therefore 
the Hawking-Penrose theorems [1^ are satisfied for this sort of singularity. However, as we evinced, the slow-roll 
parameters might become singular at the singularity points, for specific values of the free parameters that govern the 
cosmological evolution. The singularity of the slow-roll parameters might indicate some sort of instability of the dy¬ 
namical system that the cosmological evolution equations constitute. This singularity of the slow-roll parameters was 
also observed in [IMl when evolution is governed by scalars, so a repeating pattern seems to underlie this feature, 
which by no chance is accidental. This issue needs to be further scrutinized, and in general all the consequences of 
finite time singularities should be fully understood in the context of classical cosmology. This is owing to the fact 
that singularities, and especially the crushing singularities, appearing in classical theories are usually indicators of 
an underlying theoretical framework yet to be found. Therefore, the “mild” singularities might be some classical 
link between the quantum and the classical theory of gravity (apart from Type IV singularity where quantum effects 
should not play an essential role). 

An interesting possibility that we did not address in this letter is to study a Type IV singular evolution caused 
by two distinct phenomenological EoS dark fluids, or even couple one fluid with dark matter. This study however 
exceeds the introductory purposes of this paper and is deferred to a future work. 

Finally, a comment is in order: The generalized models we used are not by any means just mathematical construc¬ 
tions, since such EoS models are typical for acceleration of the Universe, since their EoS parameter is around the 
value —I. Their specific form can be imposed in principle, by a specific evolution of the Universe, therefore these EoS 
are physically motivated by the early/current status of the evolution of the Universe. 
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Appendix 


Here we provide the detailed functional form of the parameters Ji and J 2 that appear in the expressions giving the 
observational indices, namely in Fq. dSID, as functions of the effective energy density, and the e—folding number. The 
detailed form of these is. 


, ^ 1 


fipm - 


5 f{p{N))f{p{N)) ( f{p{N)) y 1 p'{N) 


p{N) 


V J ^fipiN)) 




. _45/(p(iV)) 


p{N) 




V piN) ; 




- 9 fip{N)) - 


+ 3 f4/'(p(A))-7^^^^-b2 ) (/'(p(7V))- 


U'{p{N)) \ 
2 p{N) i 


-45/'(p(iV))+9 


fjpm 

p{N) 


PiN) 

{f'ipiN))y + f{p{N))f{p{N)) 


1 fipm] .(fipm] ^ 

2 p{N) p{N) ) p{N) 

2 ~ 


-2 


fiPiN))f{piN)) ( fjpm ]' 


PiN) 


( fjpm ] i 3 ( fjpjN)) ] fp'jN) 

[ p{N) ) I 2 p[N) ) \ p{N) 

n fip{N))f{p{N)) 5 ( fipjN)) ]' 

2 p{N) 2 I p{N) I 


\ PiN) / 


2fipmf"ipiN)) 


f p''iN) \ 
\ PiN) J 


\fiPiN))ripm-2l^JM 


fp'iN)Y 

\piN)J 


(3f{p{N))f{piN)) + f{p{N))r{p{N))) p{N) - 3 (/'(^(iV)))' 

-mpmfipm +_ 3 (Mm] 


PiN) 


\ PiN) J 


(59) 
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